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ABSTRACT 

A summary of the Weibull distribution and three repre- 
sentative sampling procedures to determine point and con- 
fidence interval estimates of the parameters that occur 
in the functional form of the Weibull distribution are 
presented. Following this, a model, assuming the Weibull 
distribution, is proposed which Sead have possible ap-= 
plications in analyzing Polaris Missile System Trouble 
Failure Reports to determine point Sartienitee of the 
reliability of the Polaris Missile components. 

I am indebted to Professor P. W. Zehna for his con- 
tinued patience, encouragement, and most capable guidance 
while acting as faculty advisor. I also wish to thank 
Professor W. M. Woods for his valuable assistance as | 
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CHAPTER I 
INTRODUCTION 


This thesis is concerned with summarizing the Weibull 
distribution and providing several possible methods of 
analyzing Polaris Missile System Trouble Failure Reports 
to determine the reliability of Polaris Missile com- 
ponents when the Weibull distribution is assumed to 
characterize the failure times of the components, 

The Weibull distribution, which was first proposed 
by Waloddi Weibull [1] in 1951, is a member of a class of 
distributions which characterize wear out failure; the 
longer an item has been used the greater the probability 
of failure. Wear out failure is associated with gradual 
depletion of material, fatigue, or accumulated shocks 
and other factors. 

The key assumption underlying many of the statistical 
procedures in current use for evaluating the reliability 
of components or equipment is that the failure times 
follow the exponential distribution. However, in many 
practical situations sufficient data is not usually 
available to justify the assumption that the exponential 
failure law truly characterizes the time to failure 
distribution of the equipment or components being tested. 





The Weibull distribution which includes the exponential 
distribution as a special case is a more general assump- 
tion about the time to failure law of components or 
equipment. Kao {2} has shown that the failure age 
distribution of electron tubes is better chacterized 
by the Weibull distribution with shape parameter ap- 
proximately equal to 1.7 than by the exponential dis-~ 
tribution. Also when the exponential distribution is 
assumed to be the time to failure law, M. Zelen and 

M. C. Dannemiller [3] demonstrate that four commonly 
employed sampling procedures are not robust* with re- 
spect to Weibull alternatives. 

Chapter II summarizes the distribution theory 
associated with the Weibull distribution and defines 
terms to be used in the remainder of this thesis. 
Chapter III outlines sampling plans which can be em- 
ployed to obtain estimates of the parameters from exe 
perimental’ measurements. The various estimates to 
follow are derived under each of these plans. Chapter 
IV proposes three possible methods of using Polaris 
Missile Trouble Failure Reports for determining the | 


reliability of Polaris system components, 


# Sampling procedures which are not very sensitive to 
departures from basic assumptions are called robust. 





CHAPTER II 
SUMMARY OF THE WEIBULL DISTRIBUTION 


Let X be a random variable which may take on any 
value in the continuous scale from m to infinity where 
mis a non-negative number, The three parameter Weibull 


cumulative distribution function (¢.da.f.) 1s defined as 


~ (xm) 
(1). F(x) = l-e < for x = m, where 
0 x<m 
acs scale parameter, a > 0 


Le) 


shape parameter, b > 0 
m = location parameter”, m 20 
The probability density function is defined as the 
first derivative of F(x) with respect to x, and is given by, | 


-(x-m) 
a bel a 
(2). ax) —~ f(x) = b(x-m) a _>xz2m 
; ry 
0 x<—m 


Weibull ij 1] demonstrated that yield strength of 
Bofors steel, size distribution of fly ash, and fiber 


%# Xo, m, and X,, respectively were the notations used 
by Weibull QJ] . | 





strength of Indian cotton are distributed in good agree- 
ment with the cumulative distribution function, equation 
(1). Kao (i, J has shown that the cumulative distribution 
function, equation (1), qualifies as a failure distribu- 
tion and has successfully employed the Weibull distribu-e 
tion with location parameter, m, equal to zero to chare- 

- acterize the time to failure of electron tubes. For the 
purposes of this thesis and the sampling plans which 
follow in Chapter III, the location parameter, m, will 
be assumed to be zero. This is a reasonable assumption 
since if an item is placed in operation at time zero, it 
is exposed to the risk of failure from the time it is 
first put into use, 

For completness of presentation, the moments and: 
characteristic function of the three parameter Weibull 
distribution are listed below, but a detailed discussion 
will only be presented for the two parameter sub-family 
having m=0, Giwen f(x) as the above probability den- 
sity function, the nth moment is defined by the formula 


=-(x-m)> 


a bel a 
E(x") = x2 b(x-m) e ax 
mm 


& 


and is given by 


n | k 
n\ b n-k k 
E(x") = 2, (ice m T's + 1) 
3 





The mean and variance, denoted by A and o* respectively, 


are easily determined from the above relation to be 


Ll 
5 1 
Rema | ! (+1) 


and 


4 
a 
a b 2 1 
att (b+1) - | ie +1) | 
respectively. The function T+) is the gamma 
function which is defined for every v 70 by 


Tw) = \ xV~l o-X ax 


© 
Suppose now that T is a random variable which may 


take on any value in the continuous time scale from zero 
to infinity, and the random variable T popeeeenuemt ne 

total time to failure of an item. Let us further assume 
that the two parameter Weibull probability density funce- 
tion, f(t), specifies the probability law of the random 


variable fT where y. 
(3). f(t) = (btP-l/a) exp (-t?/a), t 20; a,b, w0 


= 0 elsewhere 
Accordingly, if an item is placed in operation at time 
zero, the probability that the total life of the item 


is less than or equal to t is expressed as 





t 
P(T Ss t) = \. (btP-1/a) exp (-t?/a)at = 1-exp(-t/a) 
Therefore the cumulative distribution function of the 


random variable T is given by 


i. Pie) = ‘1 - exp (-t°/a) 
The reliability function, denoted by R(t), of 
an item is also of particular interest since it is 
desirable to know the probability that an item or 
component does not fail in the time interval (0,t). 
Thus the reliability function is the probability 
that the time to failure of an item is at least t 
time- units and is expressed by 
pe! SO / 
R(t) = P(T >t) = i. (pt°-1/a )oxp(-t?/a)dt, for t>0 
which upon integration becomes | 
(5). R(t) = exp (-t9/a) 3 t 7 0 
Note that the reliability function 1s merely the unity 


complement of the cumulative distribution function, 


equation (lh). 





The nth moment of the random variable T is defined 
by the formula 


oo 
B(T™) =  em(ptb=2/a) exp (=-t>/a) at 
and is given by 
: . | 
(6). E(T) = =} (B +1) 
The characteristic function of the random var- 


lable T is 


"K 


Diu) = E(eitT) = J 


which upon integration becomes 


ao 
git ap t 057 e=(t/a) at 
- 3 


& 


iu) = 5° (tuys g4/b Ten 


47° 
The mean time to failure of the random variable T; 
denoted by A , is defined as the first moment about 


zero and is easily determined from equation (6) to be 


1 
b 1 
(Me Ame TB tay om 
In a similar manner, the variance denoted by O* of the 


random variable T is determined by the formula, | 





E(T°) - E(T)°, to be 


e 2 
> 3 T'é 

G= a (b #1) - (Fils) 

Another function of a time to failure probability 
distribution which provides a great deal of insight 
into the probability law is the hazard function or 
instantaneous failure rate which will be denoted by 


Z(t). The instantaneous failure rate or hazard 


function is defined by 


Z(t) . Brey 


where f(t) and R(t) are the density and reliability 
functions respectively. Substitution of f(t) and 
R(t) of equations (3) and (5) in the above express- 


ion,we have 
(8). Z(t) = vot°-1/a for t>0 


It is interesting to note that, for values of the shape 
parameter greater than one, the instantaneous failure 
rate is an increasing function of time, therefore indi- 
cating that wear out occurs during the life of the item. 
But in the special case where the shape parameter equals 
ohe and the Weibull distribution simplifies to the ex- 
ponential distribution, it is observed that the instan- 





taneous failure rate is a constant equal to l/a, inde- 
pendent of age, 

The effect of varying the shape and scale para- 
meters can most readily be illustrated by graphs of 
the probability density function because the probe 
ability that the random variable T will lie in the 
interval (t, t+h) can be interpreted geometrically 
as the area under the probability density function in 
the interval from t to t+h. Figure one illustrates 
the probability density function for several values of 
the shape parameter, b, and with the scale parameter, 
a, equal to one. The effect of changing the value of 
the scale parameter is to merely squeeze or broaden 
the graph of the probability density function, i.e., 
an increase in the value of the scale parameter 
would broaden while a decrease in the value of the 
scale parameter would squeeze the graph of the probe 
ability density function together. The first derivae 
tive of f(t) with respect to t is | 


tO 
a 
f'(t) = b tO ¢ |= to (o-1) ] 
& a 


which shows that f(t) has a maximum at 
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if 6 is greater than one. But if the shape parameter 
has a value greater than zero but less than or equal 
to one, the probability density function is a monoe 
tonically decreasing function of t. Thus increase 
ing values of the shape parameter squeezes the proe- 
bability density function together. This can be 
interpreted intuitively to mean that the faster an 
item wears out, the greater the value of the shape 
parameter, b. Figure 2 illustrates the reliability . 
function, R(t), for several values of the shape 
parameter and with the scale parameter, By equal 

to one. Notice for fixed values of the scale para- 
meter, all the curves intersect at t equal to one, 
regardless of the value of the shape parameter, 

Thus far, consideration has only been given to a 
single item. It would be desirable to be able to cone 
sider a component or complex item consisting of M 
individuals or single items. Instead of assuming that 
the items are independent, it will be assumed that the 
M items of the complex item are quasi-independent in 
the following sense. If the life of the ith item is 
@ random variable T, with corresponding density func- 
tion F, (t) then let T be a random variable whose 


cumulative distribution function, G(t), is determined 


Li 
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by the equation 
1-G(t) = P(T) 7 t)P(fo 7 t) . . . P(Ty > t) 


for t greater than or equal to zero Cs] e Solving 
the above equation for G(t) yields 

mM 

a a t>4 
7 oe 
(9). G(t) = l-e » t20 
The first derivative of G(t) with respect to t de-~ 
noted by g(t) gives 
M 
@ NT Py 
— a! By 
(10). aG(t)g(t)= bat [97 ae 
dt 


azxJ/. Al 

Equations (9) and (10) above will define the pro- 
bability density and cumulative distribution func- 
tion respectively of a complex item rs] p9 LeOey A 
component consisting of more than one item which are 
quasi-independent. The reliability Sunettion of a 
complex ame ify determined by the probability state-~ 
ment, P(T>t), to be | 


‘é t 
_ Wy a4 
(11). R(t) = g(t) dt=e ter 
: x 
and the hazard function is determined from the re- 


lation Z(t). g(t) to be 
S-SCE) 


13 





” 
ti27.” ze) = )"b,e1 "1/0, 
= 


The formula for computing the nth moment of a come 


plex item is then given by 


oo , . 
ad | 
E(T™) = Ve : (byt°int a, ) exp (- ee 
o GT peal 


which can be integrated by parts to yield 


@ 


“ | a 
(13). E(T™) = \ nen"? exp (- >! tOi/a, at 
. 3 42 


Again the mean and variance will be denoted by A and cs 


respectively and are determined by equation (13) to be 


=< mM 
(44). A= ) exp (- 2- t’ t/a, )at 


a ~ ae a 
(15). 0 = Set exp (- = t t/a, )at - 2X 
: © 


al 
Equations (14) and (15) can be evaluated by numerical 
methods. 

If the Weibull family of probability distributions 
is assumed to characterize the time to failure of an 
item or component consisting of more than one item 
where the various members of the family differ only 
by the values of the parameters occuring in the func- 
tional form of the distribution, it is clear that the 


reliability, hazard, density, and cumulative distri- 


Ly 





bution functions, along with the various moments are 
functions of the parameters, a and b. Consequently, 
useful probability statements about these items can 
only be made to the extent that estimates of these 
parameters may be given. At the present time, point 
and confidence interval estimation of parameters are 
most frequently employed. 

The maximum likelihood estimates of the para- 
meters, a and b, will be derived for each of the 
sampling procedures presented in chapters III and 
IV, although confidence interval estimates will only 
be enumerated for the scale parameter, a, when the 
shape parameter, b, is assumed to be known in one of 
the sampling plans. A thorough search of recent teche- 
nical publications indicates that confidence interval 
estimates, when both parameters are unknown, have not 
been derived and the derivation of these confidence 


interval estimates is beyond the scope of this thesis. 
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CHAPTER III 
SAMPLING PLANS 


In this chapter, several methods of obtaining 
experimental data to determine the maximum 141kelihood 
estimates (M.L.E.) of the parameters that appear in the 
functional form of the two parameter Weibull probability 
distribution will be discussed. Experimental data can 
be obtained in a variety of ways, and each particular 
method of obtaining experimental measurements will be 
defined as a sampling plan. For each of the sampling 
plans presented, it is assumed that the random var- 
iable T characterizes the time to failure of an item 
‘and the density function, f(t), defined by equation 
(3), specifies the probability law of T. It is fure 
ther assumed that it is possible to make repeated 
measurements, denoted by Ty» Tos « « « » By, Of the 
random variable T, and the measurements, Tis 0 © © sty, 
are independent identically distributed random vare- 
lables, each with the same distribution as the ran- 
dom variable T. Accordingly, after selecting a ran- 
dom sample and performing an experiment, the observed 


values of the random variable T are designated by 


ty» Cos e ¢@ e ptre 


16 





In the case of a complex item composed of M ° 
individual items, a random sample of size n, will 
be taken from the item of the ith type, for 
ZZ =1, 2, « « « »Me Again the random variables 
Thy tor fixed 1/arid Farts 2,"o% . ; Ny, are 
independent identically distributed random var- 
jables. Also after performing an experiment, the 
observed values of the random variable T; will be 
denoted by tay» Cid,y © eco » Cans. Each of the 
sampling plans presented are based upon random 
samples as defined in the above paragraph and are 
restatements of the sampling plans contained in 
reference fs] ° 

The method of maximum likelihood is utilized 
to obtain estimates of the parameters, because thier 
procedure generally tends to have desirable proper- | 
ties not always shared by other methods. In brief, 
the method of maximum 14K614hood estimation consists 
of determining the values of the parameters that max- 
imize the likelihood function, denoted by L. The 
likelihood function, L, is a function of the para- 
meters in the joint density with fixed observations. 





This maximizing process is generally accompolished by 
determining the values of the parameters that maximize 
the natural logarithm of L, since ln L is an increas- 
ing function of L. The functional form of L will vary 
in general with each sampling plan. Lloyd and Lipow 
[6] and Mood Gal provide a detailed discussion of 
the maximum likelihood method. 

The maximum likelihood estimates of the para- 
meters, a and b, will be denoted by a and > res= 
pectively in each of the sampling plans discussed 
in this paragraph. These sampling plans are out- 
lined as follows: 

SAMPLING PLAN I = Fixed sample size. 

Single item case - A random sample of size n is selec- 
- and each member of the sample is tested until it 
fails. The total time to failure of each member of 
the sample is observed and recorded. If the obser- 
ved failure times of the random variable T are 

ty, tos « « « » ty» then the likelihood function 


is given by 
n s 
= fi ttea) 
as 


Upon substituting equation (3) for f(t, ) in the 
above relation and taking the natural logarithm 


18 





of L, we have 


In L=ninbd-nina x" ar 2 +,%/e 


Att A084 


Therefore, if both the parameters, a and bd, are un- 


known, the M.L.E. of these parameters are given by 


- nin 
(16). as ? s/n 


A® 6 
A n Zz 
(17). db = na/( q, £4? Int, -& Zin t,) 
respectively. Lloyd and Lipow [6] outline a pro- 
cess of iteration which will solve equations (16) 
a ; 
and (17) for a and bd. In the event that the para- 
_ meter, b, 1s assumed to be known, the M.L.E. of the 


parameter, a, under this sampling plan is given by 


(18). a = >t, n 


wie Att 
in this situation, it is also observed that a cone 
fidence interval estimate for the scale parameter 


can easily be determined using the statistic, 


_." 
(2 Ze ty )/a 
Az 4 
which has a Chi-square distribution with 2n degrees | 


of freedom, Complete details of confidence inter- 


val estimates can be found in Mood L7I ° 


? 





Complex item case = Assuming that the complex item is 
composed of M individual items, the procedure is to 
select a random sample of size Ny from the item of 
the ith type, for i= 1, 2, « « « » Me Each member 
of the ith random sample is tested until it fails, 
and the total time to failure is observed and re- 
corded for each member of the ith sample. If the 
observed failure times of the random variable Ty, 

are tay, toy. ¢ © « » tangy, then the likelihood 
function for the item of the ith type is given by 


NG 


Per t= 1, 2, « © « » Me Therefore, upon subse 
tituting equation (3) for f(t, 5) in the above re- 


lation, we obtain 


Ly = TT (oy tas )/as exp (-t44 Jaz) 
¥ 
Accordingly, the M.L.E. of the scale and shape parae- 
meters for the item of the ith type are determined, 


in the same manner as for a single item, to be 


ni A 

A b 
ign «,= 2 t,, -/n 
a> G15 Ms 
and A 
a aA A Ns db 
(20). b, =n,a,/( 243 Lin & 

na ’ 

-a ln t 

i Z/ 


ij 


13) 


20 





mapsoulvely, ror T= 1, 2, ...+, M Again, a 
process of iteration can be employed to solve for 
the mamimum likelihood estimates. In the situation 
where all of the shape parameters are assumed to be 
known, the M.L.E. of &, under this sampling plan is 
given by 


ay = z tay Pon. 


for i=l, ay o «© © » Me | 
SAMPLING PLAN II = Fixed sample size with item 
truncation. 
Single item case - A random sample of size n is 
placed on test simultaneously, and the life test 
is terminated when r members of the sample fail. 
fhe number of failures r is fixed in advance of 
the experiment, and r is an integer greater than — 
zero but less than or equal ton. Using this 
sampling plan, the first r failures times of the 
random variable T are observed and recorded as 
the failures occur. Then, if O< t, 4t,5 oe oe sg 
<t,< oo are the first r failures, Halperin (8 ] 
indicates that the likelihood function is given by 
(21). Lent/(n-r)3 Tt f(t, ) 1-7(¢,) | a” 


Az{ 





where L is merely the density function of the first 
rorder statistics. After substituting F(t) and f(t) 
of equation (3) and (4) in the above equation, the 

M.L.B. of the parameters, a and b, are determined in 


the usual manner to be given by 


raf ta | 
(22). a = Zt, + (n-r)t,, /r . 
a:/ 
“ww 


A s. 
(23). bar a/] S't,° in t,+(n-r)t,” in t 


az{ 


2S e| 


atl 


r 


Again, the values of a anes which solve equations 
(22) and (23) can be determined by an iteration 
process. iIn the event that the bnciie parameter 
is assumed to be known, the M.L.E. of the para- 
meter, a, under this sampling plan is given by 
i a 

(2h). a = ( > ty + (n-r }t_°)/r 

Ael 

Complex item case =- Assuming that the complex item 
is composed of M individual items, the procedure 
is to select a random sample of size n, from the 
item of the ith type, for 1=1, 2, . .., M The 


n, members of the ith random sample are placed on test 





simultaneously and the ith life test is terminated 
when r,; members of the ith sample fail. As before, 
the number of failures rg is fixed in advance of 
each experiment, and r, is an integer greater than 
zero but less than or equal to ny. The first ry 
failures times of the random variable T, are } 
observed and recorded as the failures occur, If 
the first ry failure times of the random variable 
Ty are given by 0 < ty, © tipo Se. . Stip, < OF ; 
then the likelihood function for the item of the 
ith type is 


Rite 
Ly = n,!/(n, - r,)! rs £(t,5) 


[ ) | ie he 
l-F(t 
ins 
for 1 = 1, 2, « e « »M- Froceeding in exactly the 
same manner as with a single item, the M.L.E. of the 


‘ 


parameters, a4 and Das of the item of the ith type are 


a he A . a 
(25). a4 = | Z, ti yt + (ny-ry)tay, a | /v, 


and 


A Na an 
(26). PE = | r,a,/ | z, bay 4 in bag a 


~*~ lice 
cE 2 b4 ras 
(n, r,) bin, In, $ip,784 i in ts, | 


23 





As before, equations (25) and (26) can be solved by an 
~~ “Aw 

iteration process for &; and Dy; for-i = 1, 2,5. 2. &@ 

If all of the shape parameters are assumed to be known, 


the M.L.E. of the parameter As is 


he 
nn 
° —- 0 dD j 


for i = 1, 2,5 6 « © Me 

SAMPLING PLAN. III - Fixed sample size with time truncation, 
Single item case - A random sample of size n is selected 
and each member of the sample is tested for a specified 
time, say t,. The time t,» where t, is greater than 

zero but less than infinity, is defined as the truncation 
time. The number of failures denoted by R that occur in 
the time interval, (0,t,),are observed and the R failure 
times of the random variable T recorded. Using this 
sampling procedure, R is of course a RendoltWiriabLe 
which may take the value r, where r = 0,1,2, . . «. ne 

If the eueerved failure times are ty »to» o 0 0 aptyy 

then the likelihood function is ~ 


n n-r 
L= (1) 7 “oy ] a-r(e,)| for t, <t). 


The likelihood function for this sampling plan is 
formally derived in Appendix A. Again, substituting 





f 


F(t) and f(t) of equations (3) and (lL) in the above 
relation, the M.L.E. of a and b are derived in the 


usual manner and determined to be given by 


‘5 A ea 
(28). a@ = | z t,> +4 (ner)t,” |v 
nN nr ~ 
- (29). b> = «= a/ | 2 ,° ln by +. 


“w me Nn 
(n-r)t,° int, -a A ts | 
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respectively. As before, an iteration process can be 
re A 

employed to solve equations (28) and (29) for a and b. 

Also, under this sampling plan the M.L.E. of the scale 


parameter is given by 


a | 
(30). 4 = |Z a + (n-r) > | /v 


when the shape parameter is assumed to be known. 
Complex item case - If the complex item ts assumed to 
be composed of M individual items, the experimental 
procedure is to select a random sample of size ny from 
the item of the ith type, rot i=1, 2, « « « sMe The - 
sample of the ith type is tested for a specified 
truncation time oe yand the number of failures agp ogee 


by = that occur in the time interval, (O, ta. ), are ob- 


f~ 





served and the R; failure times of the random variable 
T, recorded. As before, the number of failures, Ry» 
1s a random variable which may take the value Pas 


where Py = pi 2, eo e# @ #§ +e If bai Sa e @ BPE 


are the observed values of the r, failure times; then 


using the likelihood function, 


ase? D 
(ni) af £(t4 5) la-rces,) | toa > 


the M.L.E. of a, and b, are determined to be given by 
A SS ae ¥ | 
~(3T). a4 = LZ ty i t(ng-ry)ty + [v4 


~ 


—_ N y 5 big 
ADR) snpeet eerste LE: Sty “arses 


4 S$ 
+ (ny-r, ta by ln ba, we ln tas | 


for the item of the ith type, where i= 1, 2, ... M. 
If all the shape parameters are assumed to be known, 


the M.L.E. of the me ae: a, is 


: (33) “rn ate ay eel + (nyory)t, > vs | [vs 


for i=1, 2, « « e » M. 





Thus, it is readily seen that considerable cale 
culations are required to solve for a ana in each 
of the above sampling plans. But by utilizing a high 
speed electronic computer, the iteration process ree 
quired to solve for a ana is readily accompolished. 
Kao [5] outlines a graphical method of obtaining a 
first approximation for a a If these first ape 
proximations of a and b are provided to the computer 
as the initial trial, the iteration method should 
quickly converge to a solution for Qa and D In the 
next chapter, the above sampling plans, assuming the 
Parameters, a and b, are unknown, will be utilized in 
the proposed models for analyzing Polaris Missile 
System Trouble Failure Reports. It will be further 
assumed that high speed electronic computers are 
available to accompolish the calculations required 


A AN 
to solve for a and bd. 


bed 
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CHAPTER IV 
APPLICATION TO POLARIS MISSILE SYSTEM 


In this chapter, three models are proposed which 
utilize Polaris Missile System (PMS) Trouble Failure 
Report (TFR) as experimental data for determining 
M.L.E. of parameters. By following the procedures 
proposed in the models below, it may be possible 
to obtain a useful point estimate of quality indices, 
such as mean time to failure and reliability, of 
certain items composing the missile system. The 
models are applicable primarily to missile items 
for which a life history is maintained. An item 
will be considered to have a life history if: 

(1). It 1s possible to distinguish each item 
by some identification system from all other items 
of the same type. 

(2). :A permanent record is maintained of the 
total operating time to failure, 

(3). Assuming that the item does not fail, it 
is possible to estimate accurately the total operating 
time accumulated by the item at some future date. 

In each of the models, the total operating time 


to failure of an item is assumed to be a random variable. 
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T that is specified by the two parameter Weibull prob- 
ability distribution with both parameters considered 
to be unknown. Accordingly, if n items of the same 
type are selected at random, this sample is assumed 
to consist of n independent identically distributed 
random variables which are characterized by the same 
probability law as T. It is further assumed that 
none of the items have accumulated operating time 
prior to their purchase by the United States Navy. 
Also, an item, which has failed and been repaired, 
is considered to be the same as a newly purchased 
item with no accumulated operating time. In each 

of the models, the operating time of the items is 
accumulated by check outs or other reasons for op- 
erating the item, and the TFR which provides the 
total operating time to failure is used as the ob- 
served value of the random variable. Naturally, it 
is assumed: that the time to failure reported in a 
TFR is accurate. 

The first model proposed is Sampling Plan I 
‘which is discussed in Chapter III. Consequently, 
for a sample of size n, the M.L.E. of a and b are 
given by aimee ea (16) and (17), respectively, 
and a point estimate for the reliability and mean 
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time to failure can be obtained by substituting the 
calculated values of a and S Mente equations (5) and 
(7), respectively. In the case of a complex item, 
the procedure is the same, except that the values 
are obtained for the M.L.E. of a, and b, for 
Lal, 2, « « « » M, with equations (19) and (20). 
Again, a point estimate of the reliability and mean 
time to failure can be obtained by substituting the © 
values of ay anda Dy into equations (11) and (1), 
respectively. Notice that equation (1) must be -. 
evaluated by numerical methods to obtain a point 
estimate of A. The limitation of this model is 
the waiting time required to observe the nth failure, 
Prior to outlining the next two models which 
are applications of Sampling Plans II and III, it is 
necessary to present additional background concerning 
the PMS. Consider an item of a given type, say X, 
which continues to accumulate operating time by 
check outs until it fails. Because the PMS is an 
expanding program, items of type X are constantly 
being added to the system and also accumulating 
operating time. Consequently, if n items of type 


Xx are selected as a sample, each of the items could 


30 





have accumulated varying amounts of operating time 
depending upon when the item was purchased by the 
Navy. The effect of the above situation is to con- 
plicate item and time truncation sampling plans as 
discussed in Chapter III. But Sampling Plans II and 
III can still be used to obtain M.L.E. of the para- 
meters,if a life history is kept for items of type 
X,as follows. . 

The second method for obtaining point estimates 
of quality indices employs Sampling Plan II. The 
procedure is to select P samples of size Nj» where 
Jal, 2, « « e »P, from items of type X. The items 
in the jth sample are selected so that each of the 
members of this sample are expected to have accumu- 
lated approximately the same amount of operating 
time at a future date. For each of the P samples, 
fix in advance the number of observed failures, do- 
noted by rj, which will terminate the jth experiment. 
If equation (21) is used as the likelihood function 
for the jth sample of size nj and if the P samples 
are assumed to be independent, the likelihood 


function, denoted by & , of the P samples is given by 


¢ 


| P A N,-2 
jz' Art 
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Upon substituting F(t) and f(t) of equations (3) and 
(4) in the above relation, the M.L.E. of the para- 
Meters, a and b, are determined in the usual manner 


to be given by 


(34.). az Ee startin /n 
- 
(35). b= Fey Se In ty 
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Therefore, a point estimate of a quality index can be 
obtained by substituting equations (3h) and (35) into 
the equation which defines the quality index desired. 
In the case of a complex item, the method is completely 
analogous except there are M items of different types. 
Again, life histories must be matapeinadt ror each type 
of item. ° | 

The third model uses Sampling Plan III. As in: 
the second model, P samples of size ny where 
j=l, 2, « « @ »F, are selected from items of type 
X so that each member of the jth sample is expected 
to have accumulated the same amount of operating 


time at a future date. A truncation time, denoted 





by rie is selected for each of the j samples and the 
experiment is considered to be terminated after the 
truncation time has éligepee Por each of the j samples 
of size Nys where j=l, 2, « « « »Pe As in Sampling 
Plan III, the number of failures, denoted by Rays which 
occur during the truncation time, Cin in the jth 
sample is a random variable. Assuming that each of the 
P samples are independent, the likelihood function, de- 
noted by x » is | 


P 
& = Nb, 
> 


where L, is the likelihood function of Sampling 


J 
Plan III. Proceeding as in Sampling Plan III, 


the M.L.E. of a and b are given by 


al 4 = ye, 


—> 


(36). 


(37). b = r j8/ $48 iat 


Z ) 


2 | tinyer Jt, ° in be - 2S a, 1] 





In the case of a complex item the procedure is com-= 
pletely analogous. Therefore, point estimates of 
quality indices can be determined in both the single 
and complex item case by substituting the equations 
for the M.L.E. of a and b into the equation for the 
desired quality index. 

In the last two models, the selection of the 
members of the jth sample should occur naturally 
gince a group of items of a given type is usually 
issued to a submarine at the same time. These 
items should accumulate operating time at the same 
rate through check outs until failure. The time 
when each of the items is issued can be determined 
by examining the life history of the items; con- 
sequently, the members of the jth sample can be 


determined for either model two or three. 
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APPENDIX A 


In _ Plan III, n items are tested until a 
specified time, to has elasped. The number of failures 
denoted by R that occur in the time interval (0,t,) are 
observed and the R failure times recorded. Using this 
sampling procedure, R is a random variable which may 
take the value, r, where r=0, 1, « « « » ne AS 
before, equations (3) and () define the probability 
‘density and cumulative distribution functions res- 
pectively of the random variable T. The discussion 
which follows is a formal derivative of the likelihood 
function for Sampling Plan III. 

Maximum likelihood estimators are functions of 
the data obtained from an experiment. In this case, 
the experiment is to specify a truncation time be 
and observe the number of failures, R, and the 
failure times, T,, which occur in the time in- 
terval (0,t,); consequently, the joint density 
function of the number of failures, R, and the 
times to failure, T,, given that the failures 
occur in the time interval (0,t,) would be ap- 
propriate for determining the likelihood function. 


36 





Then 


Let A be defined as the event, 

(T]<t), To< to, . . . , TRe tp) where the Ts 

are not ordered and ty< Bour0r l=), 2. cee 

Let B be defined as the event, (R=r), where 

Re Ope ly « oe 5 Be 

Let C be defined as the event that the failures 

which occur do so before time toe 

Let the event, AB, be expressed as (T,< CorecesTn <t,) 
the problem is to obtain the distribution ws 


which represents the probability statement P(ABIC). 


Using conditional probability relations, it is ob- 


serve 


where 


ad that 


P(ABIC) = E( ABC) : {Bo} = P(A\BC)P(BIc) 


P(T,< ty» eo @ @ yT p< tpiBlc) 


P(A\BC) 





a) P| (P< ti, eee sTyK<ty, Pry 7 to geee yin7to i) Bf | 





P(BOC) 


" ie tia) fr-Pce } tes a 
(r\Tt . 4 for t,<t, 
(ee) * Pro 
_ [ee 
for ty<t, 
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frag * 
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and 


p(sle) = (f\Regy | [-F(to)] O 


therefore 


At 
F(tastorcecsty»Re r)= (AYN P(t) J 1-F(to) | n-r 
for ty<to 


Accordingly, the joint density function is obtained 
by taking the partial derivatives of the above dis- 
tribution function with respect to ts for 1=1,...,r 


and is given formally by 


Eltysecc,stmpsR® r) = (Vie (ey) 1-F (to) } n-r 
Ax | 


for ts a 


38 





























